SINGULARITIES OF THE SUSCEPTIBILITY OF AN SRB MEASURE 
IN THE PRESENCE OF STABLE-UNSTABLE TANGENCIES. * 



by David Ruellef. 



Abstract. Let p be an SRB (or "physical"), measure for the 
discrete time evolution given by a map f, and let p(A) denote 
the expectation value of a smooth function A. If f depends on 
a parameter, the derivative Sp(A) of p(A) with respect to the 
parameter is formally given by the value of the so-called suscep- 
tibility function ^(z) at z = 1. When f is a uniformly hyperbolic 
diffeomorphism, it has been proved that the power series *B(z) 
has a radius of convergence r(\I/) > 1, and that Sp(A) = ^(1), 
but it is known that r(^>) < 1 in some other cases. One reason 
why f may fail to be uniformly hyperbolic is if there are tangen- 
cies between the stable and unstable manifolds for (f,p). The 
present paper gives a crude, nonrigorous, analysis of this situ- 
ation in terms of the Hausdorff dimension d of p in the stable 
direction. We find that the tangencies produce singularities of 
for \z\ <lifd< 1/2, but only for \z\ > 1 if d > 1/2. In 
particular, if d > 1/2 we may hope that ^(1) makes sense, and 
the derivative Sp(A) = ^(1) lias thus a chance to be defined. 
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0. Introduction. 



Let / be a diffeomorphism of the compact manifold M, and p an SRB measure 1 for 
/. The derivative dxp(A) of the map / H> p in the direction of the smooth vector field 2 
X, evaluated at the smooth real function A, can be formally computed to be the value at 



We shall call \& the susceptibility 3 . 

In the uniformly hyperbolic case (i.e., if the support of p is a mixing Axiom A attractor 
for /), \1/ has a radius of convergence r(\l/) > 1. One can furthermore prove that / i— >■ p 
is differentiate and that its derivative is given by \1/(1) 4 . In nonuniformly hyperbolic 
situations these assertions may fail: r(\l/) may be < 1, and / (->■ p is presumed to be 
nondifferentiable 5 . 

The above results suggest two problems: I. proving that r(\l/) < 1, > 1, or > 1 in cases 
of some generality, and II. relating the derivative of / i— >■ p to ^(1) when this quantity is 
defined. The present note is about the first problem, and presents a nonrigorous study of 
the singularities of ^ which may occur as a result of tangencies, i.e., tangencies of stable 
and unstable manifolds for the system (/, p), assumed to have no zero Lyapunov exponent. 
(The existence of tangencies excludes uniform hyperbolicity). Our study is not rigorous, 
but suggests that r(\I/) < 1 if the partial Hausdorff dimension d of p in the stable direction 
is < |, while r(tf) > 1 if d > 1/2. This opens the possibility that, for some fat tangencies 
(d sufficiently large), ^(1) is well defined. In that case, a derivative of / \-> p may exist, 
with applications to the physical theory of linear response 6 . 

1 For a discussion of SRB measures (Sinai-Ruelle-Bowen) see for instance [10], [27], 
[2], and references given there. For recent work analyzing SRB measures for a class of 
noninvertible maps, see [1]. 

2 If we replace / : x i— > fx by x i— > fx + X(fx), then p is replaced by p + SxP to first 
order in X. The derivative of / i— > p in the direction of X , evaluated at A, is dxp(A). 

3 The physical susceptibility is defined for a continuous time dynamical system, and 
is a function of the frequency u. For the discrete time dynamics considered here, the 
susceptibility would bewi-> ^>(e luJ ) : but for simplicity we call ^ the susceptibility. 

4 The differentiability of / t— > p has been established in [9], the inequality r(^) > 1 
and the identity 5xp(A) = \I/ ( 1) are proved in [16]. There are corresponding results for 
hyperbolic flows [18], [3], and generalizations to partially hyperbolic systems [5]. 

5 r(\l/) < 1 has been proved for certain (noninvertible) unimodal maps of the interval 
[17], [8], see also [19] and work in progress by Baladi and Smania. The analysis in [19] 
strongly suggests that for a certain class of unimodal maps, the function / i-> p is nondif- 
ferentiable, even in the (weak) Whitney sense. There is also numerical evidence [4] that 
r(\l/) < 1 for the classical Henon attractor. For recent work on Henon-like diffeomorphisms, 
see [14]. 

6 A basic physical article on linear response is [21]. A review of linear response for 
dynamical systems is given in [20]. 
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1. Example: volume preserving diffeomorphisms. 

Let £ be equivalent to Lebesgue measure on M, and let the /-invariant probability 
measure p be the restriction of £ to a certain open set S C M [similarly, we may also 
consider the situation where / acts on R m , and S is a bounded open set in R m ] . If either 
suppA" C 5, or suppA C 5, we may write 

oo „ oo „ 

= $> n / i( - dx ^ X ( x ) -d x (Aof n ) = -J2 zH ^ dx ) [di^X(x)]A(rx) 

n=0 ^ S n=0 ^ S 

and therefore r($) > 1. 

If we furthermore suppose that either suppA C 5", or suppA C S and £{A) = 0, we 

have 

(/, p) exponentially mixing =>- r(^>) > 1 

Volume preserving Anosov diffeomorphisms satisfy this condition, and the same is true of 
the time 1 map of an exponentially mixing volume preserving Anosov flow (which is not 
uniformly hyperbolic). Can exponential mixing happen for non- Anosov area preserving 
diffeomorphisms in 2 dimensions? We shall now see that mixing already implies that *(1) 
is well defined, and 5xp(A) = *(1) when the derivative 5x is taken along diffeomorphisms 
preserving a (parameter dependent) volume. 

For simplicity we discuss the case S = M. Let p be a probability measure equivalent 
to Lebesgue measure on the compact manifold M. Denote by p the density of p with 
respect to Lebesgue measure on some charts. Thus (f*p)~ = p o f~ x j 'J o f~ x where 
J(x) = \det(D x f)\. Suppose now that /, p depend smoothly on a parameter, and denote 
the derivative with respect to the parameter by a prime. In particular f = Io/, J'(x) = 
J(x)[divX(fx)]. 

Writing p\ = f*p we have p\ = (p/ J) o or p(x) = J{x){p\(fx)), hence 
p'(x) = J(x)[p[(fx) + d fx ~ Pl ■ X(fx)) + J(x)[divX(/x)]p!(/x) 



or 



(p'/j) o /- 1 = p[ + dp 1 -x + [divx]p! = p[ + div(p;x) = p[ + [div^xjp! 

hence 

J dxp'(x)A(fx) = J dx p' 1 (x)A(x) + J dx p 1 (x)[d\Vp 1 X(x)]A(x) 
Imposing the invariance condition p = f* p, we have thus 

dxp'(x)A(f N+1 x) = j dxp'(x)A(x) + J2 [ dxp(x)[dwpX(x)]A( f n x) (2) 



Note that J dx p{x) = 1 implies J p(dx) (p'/p)(x) = J dxp'(x) = 0. Therefore, imposing 
mixing gives that 



/ 



dxp'(x)A(f 



N+l 




tends to when N — > oo. Equation (2) now implies that 



(z) = -J2 zU / dxp(x)[divpX(x)]A(f n x) 



n=0 



is well defined for z = 1, and J dxp'(x)A(x) = 

Conclusion: Suppose that p is f-ergodic, with density p, and that (f, p) is mixing 
on a function space S containing p'/p and A, then ^(1) is well defined, and Sxp(A) = 
f dxp'{x)A(x) = *(1). 

2. Computer simulations. 

It is accepted that, using a computer, one can approximate numerically an SRB mea- 
sure by a time average: 



for large N\ — N (and iVo moderately large); the idea is that the computed orbit f n x 
is noisy because of roundoff errors, and that this noisy orbit has an SRB time average 7 . 
The Lyapunov exponents, and the coefficient L + introduced below, can also in principle 
be determined numerically. It is therefore possible to estimate r(\l/) in particular cases, 
and to test the relations proposed above between the stable dimension d of p and the 
convergence radius r(\l/) in the presence of tangencies. For example, let dim M = 2, and 
let the Lyapunov exponents A_, A + of (/,p) satisfy A_ < < A + , so that 8 d = A+/|A_|. 
Does the presence of tangencies together with A+/|A_| < 1/2 imply r(^>) < 1? (This 
appears to be the case for the classical Henon attractor). Does A+/|A_| > 1/2 imply 
r(^) > 1? Are there examples with tangencies and r(^) > 1? 

3. Singularities of \1/ in the presence of tangencies. 

It is readily seen that the power series (1) defining the susceptibility has a radius of 
convergence r(^) > 0. Tangencies between stable and unstable manifolds for (/, p) are 
expected to produce singularities of thus limiting r(^>). A difficulty of the problem 
is that the set of points of tangency has measure zero. Note in this respect that the 
angle between stable and unstable manifolds is defined p-a.e., and that the a.e. range of 
this angle determines if tangencies are allowed or not. A similar comment can be made 
for higher order contacts of the stable- unstable manifolds. It is reasonable to exclude 

7 See [15] and, for example, Eckmann and Ruelle [6]. 
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those higher order contacts which (given the dimension of M) are nongeneric if the stable 
and unstable manifolds are regarded as independent. At a generic tangency point O, the 
unstable manifold is folded in a way which is basically 2-dimensional (corresponding to 
variables x,y introduced below). Along the orbit (f n O) we have folds which are sharper 
and sharper as n — > oo. This exponential sharpening of the folds, combined with the 
derivative d x in (1), produces the singularities of with \z\ < 1 which we want to 

study. 

One can prove that r(^) < 1 for certain unimodal maps of the interval 9 ; these maps 
are non-invertible and give a degenerate example of tangencies that is relatively accessible 
to mathematical study. In what follows we discuss a crude imitation of the 1-dimensional 
situation for higher-dimensional diffeomorphisms. In the case of unimodal maps of the 
interval with an ergodic measure p absolutely continuous with respect to Lebesgue, the 
density of p has spikes ~ \x — / n c| -1 / 2 on one side of the points f n c of the postcritical 
orbit. These spikes are at the origin of the singularities of inside of the unit circle. 
Instead of an individual postcritical point, we find for higher dimensional diffeomorphisms 
a family of tangencies of stable and unstable manifolds: think of a pile of (local) unstable 
manifolds (with tangencies) carrying part of the measure p. Morally, this means that 
the spikes are "spread out" or "smoothed" (corresponding to integration over a measure 
transverse to the unstable manifolds). This smoothing may give weaker singularities of ^ 
(i.e., larger r(\I/)). 

Let us choose local coordinates (x, I,(/,F)gRx R s_1 x R x R u_1 such that the s- 
dimensional stable manifolds are given by (y, Y) = const., and the local unstable manifold 
U through O is given by x = ay 2 , X = 0, Y = (for definiteness we take a > 0). The con- 
ditional measure of p on U is thus ~ A(dx dX dy dY) = 8{x — ay 2 )S(X)S(Y)dx dX dy dY. 
One can argue that the variable Y does not play an important role in the present discus- 
sion, and we shall omit it, which amounts to taking u = 1. Using similar local coordinates 
near fO, we assume that the map / has the form 

{x,X,y)^(e L +x,e A X,e L -y) 

where L_ < 0, L + > 0, and e A is a contraction (stronger than that given by e L ~). 

The assumption that the unstable manifolds are parallel affine manifolds is crude, and 
so is the assumption that L_,L + , and A are constant coefficients. [One might think of 
L_, L_|_, and A as Lyapunov exponents. But L + , the only one of these coefficients to appear 
in the final formulas, is really the mean rate of expansion along a forward orbit {f n O) of 
tangencies]. These crude assumptions may be in part justified by the fact that we are 
looking for the leading singular behavior associated with a subset of unstable manifolds. 
We shall use informally the notation ps (approximately equal to) and ~ (approximately 
proportional to) in trying to find the leading singular behavior. 

The contribution of the conditional measure A of p on the piece U of unstable manifold 

is 

oo „ 

V A (z)~J2 zn / MdxdX dy)X(x,X,y) ■ d (x ,x,y)(Ao f n ) 

n=0 •* 

9 See footnote 5. 
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Singularities for \z\ < 1 can only come from the component Xi of X in the ^-direction, 
giving 

oo „ 

V A (z)~J2z n / dxdySix-ay^d^oD 

n=0 "* 

oo „ 

« ^(ze L +) n Xi(0) / (f* n 5(x-ay 2 )dxdy)A' 1 (x) 

n=0 "* 

where A[(x) is the derivative of Ai(x), which is A evaluated in the coordinates (a;, 0,0) 
centered at f n O. 

One has 

S(x - ay 2 ) = -±=[8(y - ^L) + 8(y + &\ 
l^Jax y/a 

hence 

p—nL+/2 rz rz 

r n (8(x - ay 2 )dx dy) = -^=^[S(y - + % + ^—)}dx dy 



where a n = ae nL +e 2nL - . Therefore 

- £{ Z e^*r*M r 0ff (*) 

so that r(\l/ A ) = e _L+ / 2 < 1. This result is in agreement with that obtained with the 
spikes of the invariant density for unimodal maps in 1 dimension (which are limiting cases 
of diffeomorphisms with tangencies) . 

Remember however that \& is defined with the measure p rather than A. Let thus V 
be part of the measure p, carried by a pile of unstable manifolds (with tangencies) near 
O, and write 

V{dxdXdy) = J 1 (dtdX)5(x-a(t,X)( y -b(t,X)) 2 -c^X))dxdy 

where the integration is over the variable £, and 7(d£ dX) is a transverse measure of p 
in the stable direction, and we assume a(£,X) > 0. It will turn out that we obtain the 
same estimate of r(\E ,r ) for different T's, and we expect that the contributions \& r to \& of 
different T's will add up convergently for \z\ < r(\l/ r ). [Such behavior was found for the 
contributions of different spikes in the unimodal case]. The most singular part of \l/ r is of 
the form 

oo „ 

= £V / T(dxdXdy)X 1 (x,X,y)d x (A 1 o f n ) 

oo „ 

* n / dX) S(x - a(£, X)(y - &(£, X)) 2 - c(£, X))Xi(x, X, y)d x A 1 {e nL + x)dx dy 



n=0 



= Y,(ze L+ ) n / 7(dtdX)5(x-a(Z,X)(y-b(Z,X)) 2 -c(Z,X))X 1 (x,X,y)A' 1 (e nL +x)dxdy 

n=0 J 

To define A±(x) = A(x,0,0) and A[(x) we have again used coordinates (x, 0, 0) centered 
at / n O. Note that A[(e nL +x) = A[([f n (x, 0, 0)]i) where [-]i denotes the first component. 
Therefore, when n is large, the functions x i-> [/ n (x, 0, 0)]i, A' 1 (e nL +x) oscillate many 
times, with a frequency ~ nL + . 

We may replace Xi(x,X, y) by = Xi(c(f, X), X, X)), and write 

1 



/ 6{x - a(£, - &(£, X)) 2 - c(£, X))dy = -=L 
J \/a(t. 



^aJ^T) y/x-c(£,X) 
where the right-hand side is replaced by if x < c(£,X). Then 

If we let 7(d£) be the image of the measure 7^ dX)X(£, X)/y/a(£, X) by (£, X) ^ £ = 
c(£, X), we obtain finally 

~ £(^ L+ ) n / /i(x)A;(e nL +a:)rfa; where h{x) = I ^=0 

n=0 V^-l 

4. Estimates when supp7 has zero Lebesgue measure. 

We assume now that supp7 has Lebesgue measure = 0. Given x £ supp 7, let 

C = max{( G supp 7 : f < x}. 

and let 7* (dry) be the image, restricted to 77 > 0, of 7(d£) by £ >->■ 77 = — £. We have 
then 

7* (dry) /" c/)'(r))drj 



h(x) , _ _ , . 

where 4>{rf) = J Q V 7*(dt) and 0(77) ~ ry d for small 77. If < a < 1 — d we let 

Ms) = /((* - O + ^"V (77)^77 = (1 - a) - £*) + v) a ~ 2 <Kv)dv 

where we have put an upper limit +00 to the integral because it does not need a cutoff. 
Therefore 

h,(x)~(l- a) ( x -t*) 1+d f°° m dt _ c(x )d+Q _! 
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A. Assuming d < 1/2 and taking a = 1/2 we may thus conclude that h(x) ^ C(x — 
£*) d-1 / 2 , hence, if I = (£* ? £**) i s a n interval of length I of the complement of supp7, we 
may estimate 

J \h(x)\ p dx « dtt- p W 2 -V = C'l 1 -* 1 ' 2 -® 

ifl < p < (1/2 — <f) _1 . By scaling we assume that the number of intervals 7 with |/| £ 
is ~ £~ d . We have thus 

L/2-d) 



/ \h{x)\ p dx~^l 1 - d - p ( 1 / 

J o 



lip> 1, and 1— d— p(l/2— d) > 0, we have thus h G L p (and the bound p < (l—d)/(l/2—d) 
appears best possible). Let l/p + 1/q = 1, then the Fourier transform h is in L q . We have 

g 1-d ~l-d {) 

(and this bound appears best possible). If \h(s)\ ~ s~ l for large s, we need tq > 1, i.e., 
t > 1/q if 1/q satisfies (3), i.e., 

1/2 

|/i(s)|~s _t with t>T^— : 

1 — a 

We come now to the estimation of J h(x)A' 1 (e nL+ x) dx where, for large n, x i— >■ A' 1 (e nL+ x) 
is rapidly oscillating with frequency ~ nL+. Since we are interested in the most singular 
part of h, we may replace it by a function with compact support. Because A[ is a derivative, 
there is no zero-frequency contribution to the integral, and we have 

with a negligible contribution of higher harmonics. Therefore 

h(x)A' 1 (e nL +x)dx\ ~ \h(e nL +)\ ~ e~ tnL + < exp[-n^-L+] 
and the bound again appears best possible, so that ty\ (z) converges for 

\z\ < exp[-(l - = exp[ '^—j- -L+] 

i.e., r(\l/f) = exp[— (1/2 — d)L + /(l — d)], and a reasonable guess would appear to be 

r(#) =exp[-^-^-L + ] 
[hence e" L +/ 2 < r(*) < exp[(d - 1/2)L+] < 1]. 



B. Assuming 1/2 < d < 1, we write ft (which is the convolution product 7 * (•) 1//2 ) 

as 

ft~7 *(-) a_1 *(-) /3_1 
where a, /3 > 0, a + /3 = 1/2, d + a < 1, or /3 = 1/2 — a, < a < 1 — d. We have thus 

ft = h x * where fti = 7 * (O"" 1 

and we have seen that h\{x) « C a (x — "We find as in A. that we can take 

|fti(s)| ~ s~ l with £ > a/(l — d), hence 

ft( S ) ~ s -"/(l-d) s -P = s -l/2-ad/(l-d) 

so that 

1 1 ft(x)A;(e" L +x)^| ~ |ft(e- L +)| < exp[-n( T ^ + i)L+] 
and ^i(z) converges for 

\z\ < exp^ + 1 _ 1)L+] = exp[( _^L _ i )L+] 
where we may let a — > 1 — <i, hence we may estimate 

r(*ft > expfly^ - = exp[(rf - 1/2)L + ] 

which is > 1. In fact r(^Y) = exp[(d — 1/2) L+] is a reasonable guess. 

The convergence radius r(^) now depends on the behavior of (/, p) away from tangen- 
cies, and we may expect that the derivative d x in (1) plays a less important role. Therefore 
r(^) is expected to depend on the mixing properties of (/, p), over which one has some 
control [25], [26]. One may thus hope that r(*) > 1, or even r(*) > 1, and that *(1) 
is well defined. The situation where the set of tangencies is large (d > 1/2) reminds one 
of Newhouse's study of persistent tangencies (wild hyperbolic sets, infinitely many sinks, 
see [11], [12], [13]). While the situation considered by Newhouse has very discontinuous 
topology, it is not unthinkable that the particular measure p behaves differentiably in some 
sense. 

C. It is plausible that the results of A. and B. remain true without the condition that 
supp7 has zero Lebesgue measure. Furthermore, if the stable dimension d of p is > 1, one 
can write d as a sum of partial dimensions 10 , and use arguments as above. One expects 
thus that the formula r(^\) > exp[(<i — 1/2)L + ] will remain correct in that case and, as 
argued in B., we may then have r(\l/) > 1 or even r(\l/) > 1. 

If we have a continuous time dynamical system (a flow) instead of discrete time dy- 
namics (a diffeomorphism) , we expect similar results in the presence of tangencies: a 

10 See [6] Section IV. D, and references given there, in particular [10]. 



9 



susceptibility function k(u>) with singularities in the upper half w-plane if d < 1/2, no 
singularity if d > 1/2, and k(0) hopefully well defined if d > 1/2. The continuous time 
dynamical situation is that most relevant for physical applications. 

5. Physical discussion. 

In this brief physically oriented discussion we shall, for simplicity, use the language of 
discrete time dynamical systems. 

We have made above a nonrigorous analysis of how tangencies between stable and 
unstable manifolds may influence the radius of convergence r(\l/) of the susceptibility func- 
tion. We have found two different regimes depending on whether the stable dimension d 
of the SRB measure p is < 1/2 or > 1/2. 

If d < 1/2 we expect r(\l/) < 1, i.e., the tangencies cause singularities of ^(z) with 
\z\ < 1. Such singularities reflect the exponential growth of small periodic perturbations 
of the dynamics (/, p) (see [20]). Experimentally, this may be visible as resonant behavior 
when a physical system is excited by a weak periodic signal: it would be of particular 
interest to study the case of hydrodynamic turbulence. 

If d > 1/2 we expect r(\E f ) > 1 and, if d > 1/2, the value \I/ (1) may be well defined. 
Since \I/ (1) is formally related to the derivative of p with respect to /, we may hope 
that this derivative exists in some sense. This would apply to physical systems not too far 
from equilibrium (at equilibrium, p has a density, and d > 1 unless all Lyapunov exponents 
vanish) with obvious application to linear response in nonequilibrium statistical mechanics. 
For large physical systems (dimM large), when there is chaos and a density of Lyapunov 
exponents can be defined, one also expects d large by the Kaplan- Yorke formula 11 , provided 
the degrees of freedom of the large system have a sufficiently strong effective interaction. 

In view of the mathematical difficulty of analyzing dynamical systems with tangencies, 
a computer-experimental study would be desirable. The situation of choice would be that of 
2-dimensional diffeomorphisms with an SRB measure p such that the Lyapunov exponents 
A_, A+ satisfy A_ < < A+. In that case we know [24] that d = A+/|A_|, and the radius of 
convergence r(\l/) is also accessible numerically. For the classical Henon attractor we have 
d < 1/2, and it appears [4] that r(\l/) < 1. In other cases, studied by Ueda and coworkers 
[22], [23], visual inspection of the computer plot of the attractor seems to indicate a large 
d, and it would be desirable to estimate r(\E f ). 

6. Infinitesimally stable ergodic measures. 

Consider the general situation of a diffeomorphism / of the compact manifold M, and 
of an ergodic measure p for / on M. We want to study formally the stability of p under 
an infinitesimal change of /. 

We shall use a space T> of smooth functions on M, with dual T>*, and a space V of 
smooth vector fields on M. \{X e V, we write X( A) = j p(dx) X(x) -d x A, so that X e V* . 
Defining T : V* V* and Tf : V V by 

(Tt)(A) = t(A o /) , ((Tf)X)(fx) = (T x f)X(x) 
11 See [6] Section IV. C, and references given there, in particular [7]. 



10 



we find that ((Tf)X) A = TX. 

Consider p + X as an infinitesimal perturbation of p (it corresponds to replacing p by 
its image under x i— >■ x + X(x)). The measure p is mapped to itself by /, while p + X is 
mapped to p + TX. We say that p is infinitesimally stable (or attracting) if (T n X)(A) — > 
exponentially 12 with n whenever X e V, A e X>. It is plausible that an infinitesimally 
stable measure must be SRB. 

We perturb ftof = f + Xof, where I 6 V. If £ G £>*, the /-invariance of p + £, 

i.e., 

(p + 0(Ao(/ + Xo/)) = ( p + 0(A) 
is then given, to first order in X, by 

| p(cfe) + X(fx) ■ d fx A] + £(A o /) = p{A) + £(A) 

or X + T£ = f, hence T n £ - T n+1 £ = T n X, hence f - T 7V + 1 ^ = ^ =o rX. Therefore, 
if p is infinitesimally stable, we obtain p + £ which is /-invariant to first order by taking 

oo 

aA) = J2(T n X)(A) 

n=0 

and £ is unique such that (T n £)(A) — > for all A e D when n — > oo. This shows that the 
linear response X i-> £ is related to infinitesimal stability. 

The above considerations apply to the uniformly hyperbolic situation where p is an 
SRB measure on an Axiom A attractor. The purpose of the present paper has been to make 
plausible the infinitesimal stability of SRB measures in a different situation where there 
are stable-unstable tangencies. [Note that ^(z) = Xl^Lo z n {T n X)(A), so that r{^>) > 1 is 
equivalent to the infinitesimal stability condition that (T n X)(A) — > exponentially]. 
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